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We test the convergence property of the chiral perturbation theory (ChPT) using a lattice QCD 
calculation of pion mass and decay constant with two dynamical quark flavors. The lattice calcu- 
lation is performed using the overlap fermion formulation, which realizes exact chiral symmetry at 
finite lattice spacing. By comparing various expansion prescriptions, we find that the chiral expan- 
sion is well saturated at the next-to-leading order (NLO) for pions lighter than ~450 MeV. Better 
convergence behavior is found in particular for a resummed expansion parameter £, with which the 
lattice data in the pion mass region 290~750 MeV can be fitted well with the next-to-next-to-leading 
order (NNLO) formulae. We obtain the results in two-flavor QCD for the low energy constants I3 
and I4 as well as the pion decay constant, the chiral condensate, and the average up and down quark 
mass. 
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Chiral perturbation theory (ChPT) is a powerful tool 
to analyze the dynamics of low energy pions The 
expansion parameter in ChPT is the pion mass (or mo- 
mentum) divided by the typical scale of the underlying 
theory such as Quantum Chromodynamics (QCD). Good 
convergence of the chiral expansion is observed for phys- 
ical pions in the analysis including the next-to-next-to- 
leading order (NNLO) for the pion-pion scattering for 
instance. In the kaon mass region, on the other hand, the 
validity of ChPT is not obvious and in fact an important 
issue in many phenomenological applications. 

Lattice QCD calculation can, in principle, be used for 
a detailed test of the convergence property of ChPT, as 
one can freely vary the quark mass, typically in the range 
m s /5 ~ m s with m s the physical strange quark mass. 
However, such a direct test has been difficult, since the 
lattice regularization of the quark action explicitly vio- 
lates flavor and/or chiral symmetry in the conventional 
formulations, such as the Wilson and staggered fermions. 
One then has to introduce additional terms with un- 
known parameters in order to describe those violations 
in ChPT, hence the test requires precise continuum ex- 
trapolation. 

The aim of this article is to provide a direct comparison 
between the ChPT predictions and lattice QCD calcula- 
tions, using the overlap fermion formulation on the lat- 
tice, that preserves exact chiral symmetry at finite lattice 
spacing a [1, 0]. With the exact chiral symmetry, the 
use of the continuum ChPT is valid to describe the lat- 



tice data at a finite lattice spacing up to Lorentz violating 
corrections; the discretization error of 0(a 2 ) affects the 
value of the Low Energy Constants (LECs) and unknown 
Lorentz violating corrections. In order to make a cleaner 
analysis, we consider two-flavor QCD in this work, leav- 
ing the similar study in 2+1-flavor QCD, which intro- 
duces much more complications, for a future work. We 
calculate the pion mass and decay constant, for which 
the NNLO calculations are available in ChPT [j| . A pre- 
liminary report of this work is found in Q . 

Lattice simulations are performed on a L? s x L t = 
16 3 x 32 lattice at a lattice spacing a = 0.1184(03)(21) fm 
determined with an input tq = 0.49 fm, the Sommer scale 
defined for the heavy quark potential. At six different 
sea quark masses m sea , covering the pion mass region 
290 MeV < m-K < 750 MeV, we generate 10,000 trajecto- 
ries, among which the calculation of the pion correlator 
is carried out at every 20 trajectories. For further details 
of the simulation we refer [7[. 

In the calculation of the pion correlator, we computed 
in advance the lowest 50 conjugate pairs of eigenmode 
of the overlap-Dirac operator on each gauge configura- 
tion and stored on the disk. Then, by using the eigen- 
modes to construct a preconditioner, the inversion of the 
overlap-Dirac operator can be done with only w 15% of 
the CPU time of the full calculation. The low-modes are 
also used to improve the statistical accuracy by averag- 
ing their contribution to the correlators over 32 source 
points distributed in each time slice. The correlators are 



2 



calculated with a point source and a smeared source; the 
pion mass m n and decay constant /„. are obtained from 
a simultaneous fit of them. 

The pion decay constant f n is defined through 
(0\Afj,\Tr(p)) — if-xPij,, where A^ is the (continuum) iso- 
triplet axial- vector current. Instead of A^, we calcu- 
late the matrix element of pseudo-scalar density P lat on 
the lattice using the PCAC relation d^A^ = 2m^ at P lat 
with m^ at the bare quark mass. Since the combination 
TO kitpiat j g n0 -|. renorma ii zec i ; n0 renormalization factor 
is needed in the calculation of /„■. This is possible only 
when the chiral symmetry is exact. The renormaliza- 
tion factor for the quark mass m q — Z^ s {2 GeV) m^' 
is calculated non-perturbatively through the RI/MOM 
scheme, with which the renormalization condition is ap- 
plied at some off-shell momentum for propagators and 
vertex functions. Such a non-perturbative calculation 
suffers from the non-trivial quark mass dependence of 
the chiral condensate. By using the calculated low-modes 
explicitly, we are able to control the mass dependence to 
determine ^'/ MOM more reliably. In the chiral limit, we 
obtain (2 GeV) = 0.838(14) (03), where the second 
error arises from a subtraction of power divergence from 
the chiral condensate. The details of this calculation will 
be given elsewhere. 

Since our numerical simulation is done on a finite vol- 
ume lattice with m T L s ~ 2.9 for the lightest sea quark, 
the finite volume effect could be significant. We make a 
correction for the finite volume effect using the estimate 
within ChPT calculated up to 0{m\ /(4 7 r/ 7r ) 4 ) 8]. The 
size of the corrections for m 2 and /„• is about 5% for 
the lightest pion mass and exponentially suppressed for 
heavier data points. In addition, there is a correction due 
to fixing the global topological charge in our simulation 
0,|3| . This leads to a finite volume effect of 0(1/V) with 
V the physical space-time volume. The correction is cal- 
culable within ChPT [l(| EH depending on the value of 
topological susceptibility \t , which we calculated in [l2[ ■ 
At NLO, the correction for m 2 is similar in magnitude 
but opposite in sign to the ordinary finite volume effect 
at the lightest pion mass, and thus almost cancels. For / ff 
the finite volume effect due to the fixed topology starts 
at NLO and therefore is a subdominant effect. Note that 
the LECs appear in the calculation of these correction 
factors. We use their phenomenological values at the 
scale of physical (charged) pion mass m v + = 139.6 MeV: 

jpbys = _ Q 4 ± Q jphya = 4 3 ± jphys = 4 4 ± q 2j 

determined at the NNLO @ and lf lys = 2.9 ± 2.4. The 
errors in these values are reflected in the following anal- 
ysis assuming a gaussian distribution. 

After applying the finite volume corrections, we first 
analyze the numerical data for m 2 jm q and /„• using the 
ChPT formulae at NLO, 

ml/m q = 25(1 + ±xlnx) + c 3 x, (1) 
I-k = /(l - xlnx) + c 4 x, (2) 

where / is the pion decay constant in the chiral limit and 
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FIG. 1: Comparison of the chiral fits including the NLO 
terms for ml/m q (top) and (bottom). Fit curves to three 
lightest data points obtained with different choices of the ex- 
pansion parameter (a;, x, and f) are shown as a function of 
ml. 

B is related to the chiral condensate. Here the expansion 
is made in terms of x = 4Bm g /(47r/) 2 . The parameters 
C3 and C4 are related to the LECs Z3 hys and ll hys , re- 
spectively. At NLO, i.e. 0(x), these expressions are 
unchanged when one replaces the expansion parameter 
x by x = 2to 2 /(47t/) 2 or £ = 2m 2 r /(47r/ 7r ) 2 , where m w 
and / w denote those at a finite quark mass. Therefore, 
in a small enough pion mass region the three expansion 
parameters should describe the lattice data equally well. 

Three fit curves (x-fit, x-fit, and £-fit) for the three 
lightest pion mass points (m, < 450 MeV) are shown in 
Figure [1] as a function of m 2 . From the plot we observe 
that the different expansion parameters seem to describe 
the three lightest points equally well; the values of x 2 /dof 
are 0.30, 0.33 and 0.66 for x-, x- and £-fits. In each fit, 
the correlation between m 2 /m g and /„• for common sea 
quark mass is taken into account. Between the x- and 
x-fit, all of the resulting fit parameters are consistent. 
Among them, B and /, the LECs at the leading order 
ChPT, are also consistent with the £-fit. This indicates 
that the NLO formulae successfully describes the data. 

The agreement among the different expansion prescrip- 
tions is lost (with the deviation greater than 3a) when 
we extend the fit range to include the next lightest data 
point at to,,- ~ 520 MeV. We, therefore, conclude that 
for these quantities the NLO ChPT may be safely ap- 
plied only below w 450 MeV. 

Another important observation from Figure [T] is that 
only the £-fit reasonably describes the data beyond the 
fitted region. With the x- and x-fits the curvature due 
to the chiral logarithm is too strong to accommodate the 
heavier data points. In fact, values of the LECs with the 
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x- and i-fits are more sensitive to the fit range than the £- 
fit. This is because /, which is significantly smaller than 
/tt of our data, enters in the definition of the expansion 
parameter. Qualitatively, by replacing m q and / by rrv^ 
and f n the higher loop effects in ChPT are effectively 
resummed and the convergence of the chiral expansion is 
improved. 

We then extend the analysis to include the NNLO 
terms. Since we found that only the £-fit reasonably de- 
scribes the data beyond m ff — 450 MeV, we perform the 
NNLO analysis using the ^-expansion in the following. 
With other expansion parameters, the NNLO fits includ- 
ing heavier mass points are unstable. At the NNLO, the 
formulae in the ^-expansion are 0] 
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FIG. 2: NNLO chiral fits using all the data points for m^jra q 
(top) and f n (bottom). Data are normalized by the value in 
the chiral limit. Solid curves show the NNLO fit, and the 
truncated expansions at NLO are shown by dashed curves. 



In the terms of £ 2 ln£, the LECs at NLO appear: f phys = 
7r i phys +8r 2 phys -151n(2V2^/P hys /m w+ ) 2 , where /P hys = 
f30.7 MeV. We input the phenomenological estimate 
[phys = _ 32 o ± 4.3 to the fit. Since the data are not 
precise enough to discriminate between £ 2 In £ and £ 2 in 
the given region of £ (0.06^0.19), fit parameters a and 
(3 partially absorb the uncertainty in Z phys . In fact, our 
final results for the LECs is insensitive to l phys . 

In Figure [H we show the NNLO fits using all the data 
points (solid curves). In these plots m 2 /m 9 and / w are 
normalized by their values in the chiral limit. As ex- 
pected from the good convergence of the £-fit even at 
NLO, the NNLO formulae nicely describe the lattice data 
in the whole data region. Wc also draw a truncation at 
the NLO level (dashed curves) but using the same fit 
parameters. The difference between the NLO truncated 
curves and the NLO fit curves to the three lightest data 
points (Figurc[T]) is explained by the presence of the terms 
In £) and£(l-5£ln£) in © and flU), respectively. 
Since the factors (f — |£ln£) and (f — 5£ ln£) are signif- 
icantly larger than I in the data region, the resulting fit 
parameters C3 and C4 in the NNLO formulae are much 
lower than those of the NLO fits. This indicates that 
the determination of the NLO LECs is quite sensitive to 
whether the NNLO terms are included in the analysis, 
while the leading order LECs are stable. 

From Figure [2] we can explicitly observe the conver- 
gence behavior of the chiral expansion. For instance, at 
the kaon mass region ~ 500 MeV, the NLO term 
contributes at a —10% (+28%) level to m 2 /m g (/„■), and 
the correction at NNLO is about +3% (+18%). At least, 
the expansion is converging (NNLO is smaller than NLO) 
for both of these quantities, but quantitatively the con- 
vergence behavior depends significantly on the quantity 



of interest. For /„. the NNLO contribution is already 
substantial at the kaon mass region. 

From the £-fit, we extract the LECs of ChPT, i.e. 
the decay constant in the chiral limit /, chiral conden- 
sate S = Bf 2 /2, and the NLO LECs Z~ phys = -c 3 /B + 
ln(2V2vr//m 7r +) 2 and Z" Phys = c 4 // + ln(2\/27r//m 7r+ ) 2 . 
For each quantity, a comparison of the results between 
the NLO and the NNLO fits is shown in Figure [3] In 
each panel, the results with 5 and 6 lightest data points 
are plotted for the NNLO fit. The correlated fits give 
X 2 /dof = 1.94 and 1.40, respectively. For the NLO fits, 
we plot results obtained with 4, 5 and 6 points to show 
the stability of the fit. The x 2 /dof is less than 1.94. The 
results for these physical quantities are consistent within 
either the NLO or the NNLO fit. On the other hand, 
as seen for T phys most prominently, there is a significant 
disagreement between NLO and NNLO. This is due to 
the large NNLO coefficients as already discussed. 

We quote our final results from the NNLO fit 
with all data points: / = 111.7(3.5)(1.0)(±§;g) MeV, 

5^(2 GeV) = [235.7(5.0) (2.0) (I 12 7 ) MeV] 3 , lf ys = 
3.38(40)(24)(^ 31 ), and lf ys = 4. 12 (35) (30) (± 31 ). 
From the value at the neutral pion mass m n o — 
135.0 MeV, we obtain the average up and down 
quark mass mud an< 4 the pion decay constant as 
m^f (2 GeV) = 4.452(81) (38) (t 2 ° 27 ) MeV and f n = 

119.6(3.0) (1.0)(t§;o) MeV - In these results, the first error 
is statistical, where the error of the renormalization con- 
stant is included in quadrature for S 1 / 3 and m u d- The 
second error is systematic due to the truncation of the 
higher order corrections, which is estimated by an order 
counting with a coefficient of ss 5 as appeared at NNLO. 
For quantities carrying mass dimensions, the third error 
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FIG. 3: Comparison of the ChPT parameters obtained from 
the NLO fit and NNLO fit. Red circles, blue squares and green 
diamonds are corresponding to the results of the fit using 
the lightest 4, 5 and all 6 data points. Error bars indicate 
statistical error. Pluses and star represent reference points. 



than the two-loop result in two-flavor ChPT 14], / = 
121.9 ± 0.7 MeV, but taking account of the scale uncer- 
tainty, which is not shown in the plot, the agreement is 
more reasonable. For / and X MS (2 GeV), we also plot 
the lattice results from our independent simulation in 
the e-regime . We observe a good agreement with the 
NNLO fits. Comparison of the LECs l% bys and Z~P hys with 
the phcnomenological values [2J also favor the NNLO fits, 
especially for l^ hys . 

With the presently available computational power, the 
chiral extrapolation is still necessary in the lattice QCD 
calculations. The consistency test of the lattice data 
with ChPT as described in this paper is crucial for reli- 
able chiral extrapolation of any physical quantities to be 
calculated on the lattice. With a two-flavor simulation 
preserving exact chiral symmetry, we demonstrate that 
the lattice data are well described with the use of the 
resummed expansion parameter £ = 2m^/(47r/ 7r ) 2 . Ex- 
tension of the analysis to the case of partially quenched 
QCD 16] and to other physical quantites, such as the 
pion form factor [ijj is on-going. Also, simulations with 
exact chiral symmetry including dynamical strange quark 
are underway fiH ]. 



is from the ambiguity in the determination of rg. We 
estimate these errors from the difference of the results 
with our input tq = 0.49 fm and that with 0.465 fm [l3|. 
The third errors for l^ hys and l^ hys reflect an ambiguity 
of choosing the renormalization scale of ChPT (47t/ or 
47r/ w ). There are other possible sources of systematic 
errors that are not reflected in the error budget. They 
include the discretization effect, remaining finite volume 
effect and the effect of missing strange quark in the sea. 

In each panel of Figure [3l we also plot reference points 
(pluses and star) for comparison. Overall, with the 
NNLO fits, we find good agreement with those refer- 
ence values. For /, our result is significantly lower 
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